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Abstract: Within a supersymmetric (SUSY) type-I seesaw framework with flavor-blind 
universal boundary conditions, we study the consequences of requiring that the observed 
baryon asymmetry of the Universe be explained by either thermal or non-thermal lepto- 
genesis. In the former case, we find that the parameter space is very constrained. In the 
bulk and stop-coannihilation regions of mSUGRA parameter space (that are consistent 
with the measured dark matter abundance), lepton flavor- violating (LFV) processes are 
accessible at MEG and future experiments. However, the very high reheat temperature of 
the Universe needed after inflation (of about 10 12 GeV) leads to a severe gravitino problem, 
which disfavors either thermal leptogenesis or neutralino dark matter. Non-thermal lepto- 
genesis in the preheating phase from SUSY flat directions relaxes the gravitino problem by 
lowering the required reheat temperature. The baryon asymmetry can then be explained 
while preserving neutralino dark matter, and for the bulk or stop-coannihilation regions 
LFV processes should be observed in current or future experiments. 



1. Introduction 



Supersymmetry is perhaps the leading possibility for physics beyond the Standard Model. 
One of its nice features is that it contains natural candidates for the observed dark matter 
in the Universe. Within the super symmetric model with minimal particle content (MSSM), 
it is customary to assume flavor-blind boundary conditions at the Grand Unified Theory 
(GUT) scale, in which case the model is referred to as mSUGRA or constrained MSSM 
(CMSSM). If the dark matter particle is the lightest neutralino, the parameter space of 
mSUGRA is very tightly constrained by the precisely determined dark matter abundance 
in the Universe [1]. 

With exact R-parity conservation, neutrinos are massless in mSUGRA. However, there 
is now overwhelming evidence that neutrinos have mass and mix; for a review see Ref. [2]. 
The simplest explanation for small neutrino masses is perhaps the type-I seesaw mecha- 
nism [3] . It is therefore natural to extend mSUGRA to allow for a seesaw mechanism and 
thus for small neutrino masses. An mSUGRA-seesaw with 5O(10)-inspired boundary con- 
ditions was recently studied in [4], where it was found that neutrinos, with their Yukawa 
couplings contributing to the running of various parameters, such as the slepton mass ma- 
trices and the trilinear couplings, substantially modify the parameter space allowed by dark 
matter. Lepton flavor violation (LFV) was then studied within the same framework, and 
the LFV rates were shown to potentially differ from existing estimates by up to two orders 
of magnitude [5]. 

In this paper we add yet another constraint to the mSUGRA-seesaw+dark-matter 
scenario, namely that the baryon asymmetry of the Universe be explained by either thermal 
or non-thermal leptogenesis. 

It is well-known that in S'O(10)-inspired scenarios where the type-I seesaw mechanism 
provides the dominant contribution to neutrino masses, thermal leptogenesis [6] typically 
fails to explain the observed baryon asymmetry of the Universe. The reason is that the 
lightest right-handed (RH) neutrino is generally too light to generate enough asymmetry [7]. 
Including flavor effects [8, 9], the situation improves since the next-to-lightest RH neutrinos 
(not accounted for in Ref. [7]), can generate a large asymmetry [10]. Nevertheless, the 
scenario remains tightly constrained [11, 12], which perfectly suits our purpose: If thermal 
leptogenesis is successful only in a very restricted part of the parameter space, then definite 
predictions for LFV rates at a given point in the mSUGRA parameter space are possible. 

Thermal leptogenesis with hierarchical RH neutrinos requires the reheat temperature 
after inflation to be above 10 9 GeV [13-15]. In mSUGRA this poses a problem because 
of the overproduction of gravitinos [16]. This tension is partially alleviated if the grav- 
itino is heavier than 30 TeV, as the reheat temperature is then allowed to be as high as 
10 10 GeV [17]. As we show, such a reheat temperature is not high enough to allow for 
thermal leptogenesis from the next-to-lightest RH neutrino decays. The consequence is 
that thermal leptogenesis in our mSUGRA SO (lO)-inspired framework is inconsistent with 
neutralino dark matter (or even gravitino dark matter). 

The alternative possibility of non-thermal leptogenesis, either at reheating from infla- 
ton decay [18, 19], or at preheating [20] allows for lower reheat temperatures than thermal 
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leptogenesis. We employ the mechanism of instant preheating [21] from SUSY fiat direc- 
tions, as presented in [22]. We show that this mechanism is able to successfully explain the 
baryon asymmetry of the Universe, while maintaining the viability of the neutralino or the 
gravitino as dark matter candidates. Our predictions for LFV rates turn out to be close 
to the current bounds for mSUGRA points in the bulk region. 

It is worth mentioning that a mixed type-I + type-II [23] seesaw mechanism can be 
naturally obtained within 50(10), and leptogenesis becomes much easier [11, 24]. However, 
for the sake of minimalism and the predictiveness, we limit ourselves to a dominant type-I 
case only. 

In Section 2 we describe the framework in which we work. In Section 3 we review 
thermal leptogenesis, introducing all the necessary tools for our computation. We also 
show the numerical results for the predicted LFV rates, and comment on the gravitino 
problem. In Section 4 we perform the same analysis with non-thermal leptogenesis at 
preheating. In Section 5 we summarize our findings and conclude. 

2. SUSY-seesaw and SO(10) GUTs 

We consider the following superpotential for the MSSM augmented by singlet right-handed 
neutrinos Nf: 

f = /mssm + {f„) aj eabL a a H b u NZ + l -{M. N ) tj NfN^ , (2.1) 

where a is the lepton flavor index, i,j are generation indices, a,b are SU{2)l doublet 
indices, e a b is the totally antisymmetric tensor with e\2 = 1, and the superscript c denotes 
charge conjugation. Here, /mssm is the MSSM superpotential, L and H u are, respectively, 
the lepton doublet and up-Higgs superfields, and Mat is the Majorana mass matrix for 
the (heavy) right-handed neutrinos. At energy scales above Mjv, the light neutrino mass 
matrix is given by the type-I seesaw formula [3] , 

M v = -i v m N Xv 2 u , (2.2) 

where v u is the vacuum expectation value (VEV) of the neutral component of the up- 
type Higgs doublet H u . We denote the eigenvalues of the light neutrino mass matrix by 
m u i, i = 1,2,3. In the limit in which all RH neutrinos are decoupled, the light neutrino 
mass matrix is M. v = —nv^; n is the coupling matrix of the dimension-5 effective operator 
generated by RH neutrinos, which is determined by matching conditions at the RH neutrino 
decoupling thresholds. The matrix Ai v is diagonalized (in the basis where charged leptons 
are diagonal) by the MNS matrix that can be parameterized by the mixing angles #12, O23, 
and #13, the Dirac phase 5, and two Majorana phases 4>i an d 4>2 (see Ref. [5] for our 
convention) . 

Inspired by 50(10) GUTs, we introduce the vector R vu = (Ri, R2, R3), with strictly 
positive entries, which relates the diagonal up-type quark Yukawa couplings f„ iag to the 
diagonal neutrino Yukawa couplings f^ iag at the GUT scale: 

(f* ag ). =^(f u diag ). . (2.3) 
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Since we will always assume R4 < 0(5), a highly hierarchical pattern f^g ag 3> f^ ag 3> f^ ag 
is obtained. Note that the R^s are in general all different, and in the minimal 50(10) 
scenario, the range is typically 1 < Ri < 3. However, higher values Ri > 5 can be easily 
achieved with non-renormalizable operators. The condition in Eq. ( |2.3| ) corresponds to 
an extension of the "small mixing" scenario presented in [5]. We will not consider the 
"large mixing" case here, since large regions of the parameter space are already excluded 
by existing bounds on r — > /try [5] . 

Equation ( |2.3D implies that the RH neutrinos have a very strong hierarchy. To see 
this, assume tribimaximal mixing for the light neutrinos, and neglect the small CKM-type 
mixing in f u . For a normal hierarchy of light neutrinos, m v \ <C m U 2 <C rn U 3, one obtains [5]: 

M Nl ~ ^R\ , M N , 2 ~ 2 ™^Rl , M N:i ~ -^Lr* , (2.4) 

whereas for the inverted mass hierarchy {m ui ~ m U2 3> rn U3 ), 

M Nl ~ ^i?? , M^ 2 ~ Mflg , ~ ^-flg . (2.5) 
m V2 om Ul 2m U3 

From the above scaling behavior, we see that a quasi-degenerate spectrum (m ui ~ m U2 ~ 



m 



"3 



would require the lightest Majorana mass to be in the 10 2 -10 3 GeV range with signif- 
icant L-R mixing in the sneutrino sector, which we disregard because it would substantially 
complicate the sneutrino mass spectrum and phenomenology. Moreover, since the next- 
to-lightest RH neutrino is also lighter than in the case of the normal hierarchy, successful 
thermal leptogenesis is rendered more difficult. The inverse hierarchical case would require 
the heaviest Majorana mass to be of order 10 17 GeV, which is well above the GUT scale. 
This type of spectrum also suffers from instabilities under very small changes to M^r and 
RGE evolution [25]. For all these reasons we choose to focus on the normal hierarchy of 
light neutrinos, m v \ <C m V 2 <C rn v ^. 



3. Thermal leptogenesis 

Thermal leptogenesis is one of the most popular mechanisms to explain the observed baryon 
asymmetry of the Universe [6] . The crucial parameters for leptogenesis are the CP asym- 
metries £i a and the washout parameters Ki a . The CP asymmetry from the decay of the 
heavy (s)neutrino iVj (iVj) into a (s)lepton of flavor a is in full generality given by [26] 
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where Xi = M%./M% and 
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In the SUSY limit, the decay width is 



r(Ni -> £ a H u ) + r(Ni -> ^ Q fft) = r(iVi -> + r(iv 4 



r(iV* -> £ a H u ) = T(Ni -> = i^^M^. 

07T 



We can then define 



F(T = M Ni ) ~ m*M Ni l[tM ' [6 - 6) 

where m* ~ (1.56 x eV) sin 2 /3. 

The baryon asymmetry is obtained by solving a set of coupled Boltzmann equations 
as given for instance in Ref. [11]. However, we use convenient semi-analytical expressions 
for the final baryon asymmetry. The quantity that describes how efficiently the asymmetry 
is produced is the efficiency factor k, which is a function of Ki a . For an initial thermal 
abundance of RH (s)neutrinos, n is given by [15], 1 



n{Kia) 



K ia z B {2Ki 
where 



1 — exp 



2K ia z B (2K ia ) 



(3.4) 



z B (K) ~ 2 + 4K 0A3 exp f-^rj • ( 3 - 5 ) 

With a vanishing initial abundance of RH (s)neutrinos, a different result ensues. A fit valid 
both in the weak washout (Ki a < 3) and in the strong washout regime {Ki& > 3), was 
obtained in [27]: 



K(K h 



2.6 \ (K ic 



K in J V 0.06 



1.16" 
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(3.6) 



For the hierarchical mass spectrum of the RH (s)neutrinos <C M^ 2 <C Mjv 3 , 

the asymmetry production from each RH (s) neutrino can be considered separately, and 
eventually summed to obtain the final asymmetry. There are three mass ranges that need 
to be considered: the three-flavor regime for M Ni < (1 + tan 2 /3) x 10 9 GeV, the two-flavor 
(e + n and r) regime for (1 + tan 2 (3) x 10 9 GeV < M Ni < (1 + tan 2 /3) x 10 12 GeV and the 
unflavored regime for M Nt > (1 + tan 2 f3) x 10 12 GeV [8, 9]. 

It is well-known that with SO(10)-inspired mass relations such as in Eq. (|2.3|), the 
lightest RH (s)neutrino Ni (Ni) typically fails to produce enough asymmetry because 
is usually predicted to be much smaller than the Davidson-Ibarra bound [13] of about 
10 9 GeV for successful leptogenesis. Note however that tiny regions in the parameter space 
exist where RH neutrino masses are quasi-degenerate, Mn x ~ M^ 2 ~ Mjv 3 , in which case 
the CP asymmetry can be dramatically enhanced [28] and leptogenesis is possible [7]. We 
do not entertain this possibility any further. 



For the SM case, the replacement 2 Ki a — > Ki a must be made because of the fewer decay modes for 
each heavy particle. 
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We consider the crucial contribution to leptogenesis to arise from the next-to-lightest 
RH (s)neutrino, N 2 (N?) [10]. The asymmetry is typically produced in the two-flavor 
regime, but it is necessary to include the potential washout from N\ (Ni), which occurs in 
the three-flavor regime. We then have [12, 29] 



VB,2 



~ 0.96 x 10" 



i 2e K(K 2e + K 2 p) exp ^-^-Kie^j + e 2lI k(K 2e + K 2fl ) exp ^-^-K^ 



+ e 2t k(K 2t ) exp 



3vr 



-K lT 



(3.7) 



If the asymmetry from N 2 (N 2 ) is produced in the three- flavor regime, 



r\ B 2 ~ 0.96 x 10 



e 2e n{K 2e ) exp ^-^-K le ^ + e 2/1 K(K 2fl ) exp 
+ s 2t k(K 2 t) exp I — -^Kii 



3?r 



K 



(3.8) 



Since the asymmetry production from the heaviest RH (s)neutrino, N% (A3), occurs typi- 
cally in the unflavored regime, the CP asymmetry is suppressed by (Mn 12 /Mn 3 ) 2 , and we 
neglect it. 

The total asymmetry is given by the sum of the three RH neutrino contributions, of 
which just one is relevant, and thus 



VB ^ 1]B,2 



to be compared with the measured value [1] 



7jg MB = (6.2 ±0.15) x 10- 10 . 



(3.9) 



(3.10) 



As we shall see in the next subsection, the final baryon asymmetry produced through 
leptogenesis will be typically dependent on the initial abundance of RH (s)neutrinos; 
Eq. (3.4) is valid for an initial thermal abundance of N 2 (N 2 ) 



whereas Eq. ( |3.6| ) is valid 
for a vanishing one. Note that within our S'O(10)-inspired scenario, a thermal N 2 (N 2 )- 
abundance is very easily obtained if the Z' of U(1)b-l, which is naturally present if 5O(10) 
breaks to the left-right model, is heavier than N 2 (N 2 ) by two orders of magnitude (and 
Mz> ^$ 10 Ir, with Tr the reheat temperature) [30]. We then have that the interaction 
N 2 ^N 2 Z'^ efficiently brings the next-to-lightest RH (s)neutrino into equilibrium without 
interfering with the production mechanism from N 2 {N 2 ) decays. 

3.1 Results 

As stated in Section 2, we present results only for the normal hierarchy of light neutrinos, 
m u i <^ m u2 <^ m U 3. More precisely, we require that m u \ < m so \ = 0.009 eV. 

We use ISAJET-M to produce the neutrino and SUSY spectra [5]. The program 
implements RGE evolution in the MSSM with a type-I seesaw in full matrix form at the 
2-loop level. All sparticle masses are computed with complete 1-loop corrections and gauge 
and Yukawa coupling evolution include multiple sparticle threshold effects. The decoupling 



- 5 - 



of RH neutrinos is performed at multiple scales equal to their own running masses, and 
effects of f u on the MSSM parameters are included. It is known that large i u entries can 
significantly affect the RGE evolution with concomitant effects on the MSSM spectrum 
and the neutralino DM rates [4, 5,31]. 

In the neutrino sector, we use a top-down approach with f u and M^r input at Mgut 
and the neutrino mass matrix M. v obtained by RGE evolution. We fix ^(Mgut) using 
the SO( 10)-inspired relation Q2.3] ) and adjust Mat to produce a viable spectrum of light 
neutrinos. Since leptogenesis is sensitive to details of the MNS matrix, we scan over 5 
parameters: 4>\, <p2, 8, #13 and m v \. 

In Fig. |], we show the results for the baryon asymmetry generated by thermal lepto- 
genesis for m u \ = 0.005 eV and R v = (1,5,1) with <j>\, fa, 5 and #13 varied. Note that 
the figure displays weak scale values for all the parameters. From the bottom-right panel, 
we see that there are only a few points above the 2a lower bound on the observed baryon 
abundance [cf. Eq. ( p. 10 )], i.e., r/s > 5.9 x 10~ 10 . This clearly shows that the parameter 



space that yields successful leptogenesis is quite restricted. 

As noted in Ref. [12], the baryon asymmetry is essentially independent of R±, and 
only mildly dependent on .R3. (At the end of the section we discuss the role of R3 on the 
predictions for the LFV rates.) On the other hand, R2 is crucial for leptogenesis, in that 
it fixes the next-to-lightest RH neutrino mass scale [cf. Eq. (|2.4|)], which itself sets the size 



of the CP asymmetry [cf. Eq. (3^1)]. We confirm the finding of Ref. [12] that leptogenesis 
is only possible for R2 > 3, and set R2 = 5 in our calculations. 

To obtain more points with sufficiently large values of the baryon asymmetry, we now 
focus on a restricted parameter space. In Fig. g we show results for rjs in the parameter 
space, 2<</>i<4, — 1<02<1,1<5<5 and < #13 < 0.2, with m v \ varied. We obtain 
many allowed points, most of which have in common that the asymmetry is produced 
in the e flavor, and where K2 e + i^2u ~ 1/2, in which case the efficiency factor is close 
to maximal, and depends mildly on the initial conditions. With a vanishing initial RH 
neutrino abundance the efficiency factor would be lower by a factor of 2-3 for these points, 
so that the final asymmetry would be slightly lower than that observed. Note that although 
the CP asymmetry E2 T is typically the largest one (and so is K2 T ), the asymmetry in the tau 
flavor typically suffers from a large washout from K~2 T ] the washout from N\(N\), although 
set by a relatively small K\ T < 20, can also have a large impact due to its exponential 
effect [see Eq. Q]. 

From Fig. ||, we find a lower bound on the lightest neutrino mass and on #13: 

m„i > 0.004 eV, 9 13 > 0.04. (3.11) 

Our lower bound on the lightest neutrino mass is slightly more restrictive than found in 
Ref. [12], where a non-SUSY framework and a vanishing RH neutrino abundance were 
considered. On the other hand, our lower bound on #13 agrees well with that of Ref. [12] in 
the mass region below mi = 0.009 eV. Note that the Daya Bay [32] and Double Chooz [33] 
reactor experiments are sensitive to #13 for #13 > 0.05. 

As mentioned above, with a vanishing initial N2 (A r 2)-abundance it is very difficult to 
obtain a baryon asymmetry in the allowed range. At this point, it is worth commenting on 
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R vu =(1,5,1),m v1 ~0.005eV 




2 4 6 0.05 0.1 0.15 0.2 



Figure 1: Full parameter space scan for thermal leptogenesis (thermal initial ^-abundance). The 
color code is evident from the bottom-right panel. 

the size of the theoretical errors in the computation of the baryon asymmetry. First, we have 
neglected spectator processes [34], including flavor mixing in the so-called C matrix [8], 
which can reduce the final asymmetry by up to 30% [29]. We also neglected quantum 
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fK tt =(1,5, V 



■ m v , =0.007 eV 




Figure 2: Baryon asymmetry from thermal leptogenesis (thermal initial A^-abundance) vs. #13 for 
different values of the lightest neutrino mass m v \. The horizontal line marks the 2a lower bound, 
tie > 5.9 x 1CT 10 . 

statistical factors and assumed that kinetic equilibrium holds, an approximation that is very 
good in the strong washout regime (K\ a 3> 1), but which can make a 50% difference in the 
weak washout regime {K\ a <C 1) [35]. In our study, we obtained the largest asymmetries 
when K\ a ~ 1, in which case the uncertainty is less than 10%. Finally, we did not solve 
the full quantum Boltzmann equations based on the Keldysh-Schwinger non-equilibrium 
formalism; for recent related work see Ref. [36]. Note that in this formalism, thermal 
corrections are automatically taken into account. According to Ref. [37] an enhancement 
of the CP asymmetry parameter by a factor of a few is possible. 

It is difficult to estimate the cumulative effect of all these theoretical uncertainties 
on our results. Nevertheless, a 50% uncertainty in the final asymmetry seems to be a 
fair assessment and we conclude that a vanishing initial RH neutrino abundance in this 
framework cannot be excluded. 

Next, we extract from Fig. [2] the values of the RH neutrino masses M^ 2 required for 
successful thermal leptogenesis. The results are shown in Fig. ||. These clearly point to 
a very high reheat temperature, Tr ~ M^ 2 ~ 10 12 GeV, which poses a problem as we 
explain in the next subsection. 

We now turn to LFV. It has been known for a long time that the supersymmetric 
seesaw potentially leads to large LFV rates, due largely to slepton contributions at the 
loop level [38]. It is then a quantitative question to know if the points compatible with 
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Figure 3: Values of Mn 2 for which thermal leptogenesis is successful for different values of the 
lightest neutrino mass m v \. 



thermal leptogenesis in our framework lead to observable rates in future experiments. We 
show in Table [l] the current bounds and projected sensitivities for LFV. 

So far our results have been es- 
sentially independent of the region 
of mSUGRA parameter space com- 
patible with the dark matter abu- 
dance since leptogenesis occurs at 
very high energy scales. However, 
for the LFV rates, it is of crucial 
importance. We focus on the bulk 
region, which is the most optimistic 
region for the detection of LFV. This 
region is characterized by the fol- 
lowing parameters: mo = 80 GeV, 
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Table 1: Present bounds and projected sensitivities for 
LFV processes. 



m 



1/2 



170 GeV, A Cl 



-250 GeV 



and tan/3 = 10. We also checked that points in the stop-coannihilation region (mo = 
150 GeV, m\/2 = 300 GeV, Aq = —1095 GeV and tan/3 = 5) yield rates that are no more 
than a factor of two different than in the bulk region [5]. This is not surprising given the 
fact that LFV rates are maximized when sfermions and gauginos are light, and when the 
mass hierarchy is mild between them. 
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Bulk region, R yu =(1,5, 1), m v1 =0.005 eV, r| e > 5.9x10 
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Figure 4: LFV rates for points with an asymmetry above 5.9 x 10 10 within thermal leptogenesis 
(red dots) and non-thermal leptogenesis at preheating (blue stars). 



The results for LFV rates in the bulk region are presented in Fig. ^. We find that 
the allowed points lead to predictions for LFV rates below the present exclusion limits, 
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and, except for r — > fi^n, within reach of current and future experiments. In particu- 
lar, MEG should see a fi — > cy signal if leptogenesis is realized in our framework. It 
should be noted that LFV rates are very weakly dependent on the value of the lightest 
neutrino mass m u i, and we checked that the rates are essentially unchanged in the range 
0.004 eV < m vl < 0.009 eV. 

We end this section by commenting on the dependence of LFV rates on R3. Increasing 
i?3 leads to a larger M^ 3 [see Eq. (|2.4[) 1, which in turn leads to larger Yukawa couplings 
(fv)a3 hi order to keep the neutrino mass matrix fixed. We expect LFV rates to increase 
as we increase -R3 since the third generation Yukawa contributes dominantly to the rates. 
We explicitly checked that this is the case, and found the rates for R3 = 5 to be about one 
order of magnitude larger than those shown in Fig. |||. 

3.2 Reheat temperature and the gravitino problem 

In order for leptogenesis to work in our framework, we need a very high reheat temperature 
(~ 10 12 GeV) of the Universe after inflation (see Fig. ||). Therefore, the gravitino prob- 
lem here is even more severe than in the conventional scenario, where typically a reheat 
temperature of about 10 9 GeV is sufficient. 

Let us first consider the case where the gravitino is not the lightest supersymmetric 
particle (LSP). The gravitino problem arises because the thermal production of gravitinos 
is unavoidable after inflation. The gravitino yield, I3/2, defined as the gravitino number 
density divided by the entropy density, is nearly proportional to the reheat temperature [17]: 

V 3/2 ~ 2.3 x 10~ 14 x T^ 8) (l + 0.015 logT^ 8) - 0.0009 log 2 

+1.5 x HT 14 x T § ] (l - 0.037 log Tjf } +0.0009 log 2 T®) , (3.12) 

where T^' = Tr/(10 8 GeV), m 3 / 2 is the gravitino mass, and mi/ 2 is the unified gaugino 
mass at the GUT scale. Note also that increasing the gravitino mass lowers its abundance 
until the first term in Eq. ( |3. 12 ) dominates, at which point the abundance saturates. 



Gravitinos being only gravitationally coupled, for masses below 30 TeV they typi- 
cally decay during or after Big Bang Nucleosynthesis (BBN), hence spoiling the agreement 
between observations and theory. Requiring that BBN is successful leads to stringent 
constraints on the reheat temperature as a function of the gravitino mass (see [17] and 
references therein): Tr < 10 6 GeV if the gravitino mass is lower than about 10 TeV. 

For m 3 / 2 > 30 TeV, the bound relaxes to Tr, < 10 10 GeV so that that the gravitino 
decays into the LSP (assumed to be the lightest neutralino) yielding a non-thermal con- 
tribution that saturates the observed dark matter density. This upper bound is still at 
odds with thermal leptogenesis in our framework. Therefore, we conclude that thermal 
leptogenesis is somewhat incompatible with neutralino dark matter. 

Suppose we abandon neutralino dark matter. Consistent cosmology first requires that 
the gravitino be heavier than 30 TeV so that it decays to the LSP before the onset of 
BBN. In turn, the LSP may decay to a hidden sector [50]. The lightest hidden sector 
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particle X needs to be much lighter than the LSP so that its energy density is diluted 
before matter-radiation equality: 

Q x h 2 ~ 2.8 x 10 5 x Y 3/2 (y^) • (3-13) 

Using Eq. ( |3.12| ) for a very heavy gravitino S> mi/2, we find 

n x h 2 ~ 6 x 10- 5 f 19 Tr f -^U . (3.14) 

The hidden sector particle could constitute part of the dark matter, but being warm, its 
abundance cannot exceed 5% of the total dark matter abundance [51]. 

An alternative way to circumvent the gravitino problem is if the gravitino is the (visible) 
LSP, with a mass of about 200-300 GeV, and itself decays into a much lighter hidden 
sector particle long before matter-radiation equality. In this case, the NLSP must decay 
before BBN, which can be easily achieved with hidden sector dynamics [50]. Note that the 
gravitino cannot be the dark matter particle. 

Independently of 7713/2, we find that dark matter must be explained by an external 
mechansim and particle, like the axion. 

4. Non-thermal leptogenesis at preheating 

As explained in the previous section, within our mSUGRA-seesaw framework with S'O(IO)- 
inspired boundary conditions, standard thermal leptogenesis is subject to a severe gravitino 
problem. 

Non-thermal leptogenesis allows for a low reheat temperature so that neutralino dark 
matter is viable. Non-thermal RH neutrino production can be obtained either during the 
preheating stage [20], or from inflaton decays [18, 19]. 

Here we follow the approach of Ref. [22] which does not require any additional ingre- 
dient to our framework such as a large coupling of the RH neutrinos to the inflaton. The 
idea is to use the presence of flat directions in the scalar potential (for a review see [52]) to 
enable instant preheating [21], which is a very efficient way of producing very heavy states. 
We briefly review the mechanism. 

F- and D-term flat directions can be lifted by soft supersymmetry (SUSY) breaking 
terms in our vacuum, with non-renormalizable terms in the superpotential, or with finite 
density terms in the potential [53]. If the soft SUSY-breaking term has a positive sign 
and the finite density term (proportional to the Hubble rate squared H 2 ), contributes 
negatively, the field <f> along the flat direction initially acquires a large vacuum expectation 
value (VEV) denoted by 4>q. After inflation ends, the inflaton starts oscillating at the 
minimum of its potential, while the Hubble rate falls. Once H ~ m/3, where m is a 
soft SUSY-breaking mass term, the flat direction starts moving down towards the true 
minimum at <p = 0. Following [22], we assume that the condensate involves the third 
generation quark U3, and focus on the production of the up- type scalar Higgs H u relevant 
for leptogenesis. The condensate couples to the Higgs through the term ft\<p\ 2 \H u \ 2 . If the 
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condensate passes through the origin (or sufficiently close to it), Higgses will be produced 
when adiabaticity is violated [21], i.e. mH u /m 2 Hu > 1. 

The condensate continues its motion upwards after it has passed through the origin, 
and the up-type Higgs effective mass gradually increases proportionally to ft\4>\- When 
the Higgs effective mass becomes larger than the RH neutrino mass, it promptly decays to 
RH neutrinos. The heaviest particles that can be produced through this mechanism have 
a mass [22] 

/ I i I \ 1/2 / ~ \ 1/2 
M max „ 4 x 1Q 12 QeV ( ™ ) , (4.1) 



VMpiy ViOOGeV, 

which is large enough for our purposes. Accounting for the fact that reheating occurs after 
leptogenesis, a large dilution factor must be included. Assuming that all up-type Higgses 
decay into RH neutrinos, the final result for the baryon asymmetry is given by [22] 



which means that for the canonical choice of parameters the efficiency factor is about 
7 x 1(T 3 . 

For a highly hierarchical RH neutrino mass spectrum, and small , we need to again 
consider the production of asymmetry from N% . Specifically, we need the up- type Higgs to 
promptly decay into N2 rather than N%, which is possible if [22] 

M«, < ( I'' 2 ■ (4.3) 



The condition can be easily satisfied if |(fi/) Q i| "C 1. 

Finally, we need to take into account the washout by N\ inverse decays of the asym- 



metry produced by N2. Using Eq. fl4.2|) we find 



VB 



2 -7xl0- 5 ^e 2a exp f-^K la j , (4.4) 



where we have taken Tr = 10 8 GeV, m = 100 GeV and <fro = Mp\. 

Before turning to the numerical results, we comment on another non-thermal scenario 
of leptogenesis, namely at reheating from inflaton decay [18, 19]. The final asymmetry 
obtained for an inflaton of mass 10 13 GeV is [18] 

VB, * I x 10- 5 5> a Br(0 N^) Lj*^ ) , (4-5) 

where Br(^ — > NiN{) is the branching ratio of the inflaton decay channel <p — > NiNi. We 
find that even if the branching ratio to the next-to-lightest RH neutrino N2 is unity, the 
maximum efficiency factor is a factor of 2 smaller than in the preheating case. All the 
results presented in the next subsection can then be trivially extended to the reheating 
case. 
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4.1 Results 

Similarly to the previous section, we explore the parameter space for which Eq. ( f4.4| ) re- 
produces the observed baryon asymmetry; see Fig. ||. Note that the measured baryon 
asymmetry can be generated for more of the parameter space than with thermal lepto- 
genesis. The reason is simply that the asymmetry depends on fewer paremeters than in 
the thermal case. In particular, the asymmetry does not depend on Ki& as was the case 



for thermal leptogenesis [see Eqs. (44) and (13.71 )]. Therefore, non-thermal leptogenesis is 
possible even for very large values of the washout parameter K<i a , so long as K\ a remains 
small. In the non-thermal case, the SO( 10) mass relations in Eq. (|2.3| ) and light neutrino 
masses lead us to a new part of the parameter space in which the Yukawa coupling (i u ) T 2 
(and therefore K^t) is large, while K\ r remains small. In this region, the lepton asymmetry 
is produced in the t flavor because of the large value of the CP asymmetry parameter S2 T 
and can be up to two orders of magnitude larger than the observed value. Due to the 
greater freedom in the choice of parameters, it is not surprising that neither m v \ nor #13 
are bounded from below. 

For the points in Fig. ||, we show the corresponding LFV rates in Fig. ||. As for thermal 
leptogenesis, we find that large rates are predicted for the bulk and stop-coannihilation 
regions. In particular, MEG should see a positive signal if leptogenesis is the origin of the 
baryon asymmetry. Note that the rates remain essentially unchanged under variations of 
the lightest neutrino mass m v \ within the range of interest. 

Finally, let us comment on the gravitino problem in this framework. Since the reheat 
temperature after inflation is required to be of order 10 7 -10 8 GeV, it is clearly not as severe 
as in the previous section. However, if dark matter is to be explained by the standard 
neutralino LSP, we still need a fairly heavy gravitino of about 10 TeV [17]. On the other 
hand, if the gravitino is the LSP, we must ensure that the NLSP decays before BBN, which 
can be easily achieved with small R-parity violation [54] , or with decays into a light hidden 
sector [50]. 



5. Conclusions 

We studied the implications of successful leptogenesis on the mSUGRA-seesaw parame- 
ter space with dark matter comprised of neutralinos. Guided by S'O(10)-inspired mass 
relations, we were led to hierarchical Dirac mass eigenvalues in the neutrino sector. 

We found that with thermal leptogenesis, a large enough baryon asymmetry is difficult 
to realize, and obtained lower bounds on both the lepton mixing angle #13 and the lightest 
neutrino mass m v \. The LFV rates in the bulk and stop-coannihilation regions are large 
and observable at current experiments such as MEG. However, the high reheat tempera- 
ture (~ 10 12 GeV) implies a severe overproduction of gravitinos, rendering neutralino dark 
matter and thermal leptogenesis somewhat incompatible in our framework. 

To relax the tension with gravitino overproduction, so as to not abandon neutralino 
dark matter, we explored the possibility of non-thermal leptogenesis in the preheating phase 
which relies on the mechanism of instant preheating from super symmetric flat directions 
already present in our model. We found that an efficient (non-thermal) production of 
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R vu =(1, 5, 1), m vl ~0.005 eV, T R =1(f GeV 
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Figure 5: Full parameter space scan for non-thermal leptogenesis at preheating. The color code 
is evident from the bottom-right panel. 



RH neutrinos can be achieved, even for reheat temperatures as low as 10 7 GeV. Note 
that non-thermal leptogenesis from innaton decay would lead to very similar results. The 
parameter space for successful leptogenesis at preheating is less constrained than in the 
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thermal framework, and LFV rates in the bulk and stop-coannihilation regions remain large 
and observable at current and future experiments. 
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